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ABSTRACT 

The Spitzer Adaptation of the Red-sequence Cluster Survey (SpARCS) is a deep z'-band imaging 
survey covering the Spitzer SWIRE Legacy fields designed to create the first large homogeneously- 
selected sample of massive clusters at z > 1 using an infrared adaptation of the cluster red-sequence 
method. We present an overview of the northern component of the survey which has been observed 
with CFHT /MegaCam and covers 28.3 deg 2 . The southern component of the survey was observed with 
CTIO/MOSAICII, covers 13.6 deg 2 , and is summarized in a companion paper by Wilson et al. (2008). 

We also present spectroscopic confirmation of two rich cluster candidates at z ~ 1.2. Based on Nod- 
and-Shuffle spectroscopy from GMOS-N on Gemini there are 17 and 28 confirmed cluster members in 
SpARCS J163435+402151 and SpARCS J163852+403843 which have spectroscopic redshifts of 1.1798 
and 1.1963, respectively. The clusters have velocity dispersions of 490 ± 140 km/s and 650 ± 160 
km/s, respectively which imply masses (M 200 ) of (1.0 ± 0.9) x 10 14 M 0 and (2.4 ± 1.8) x 10 14 M 0 . 
Confirmation of these candidates as bona fide massive clusters demonstrates that two-filter imaging 
is an effective, yet observationally efficient, method for selecting clusters at z > 1. 

Subject headings: infrared: galaxies 


1. INTRODUCTION 

In the nearby universe there are numerous lines of ev- 
idence suggesting that environmental processes could be 
the dominant force driving the evolution of the galaxy 
population. Properties such as star formation rate (SFR, 
e.g., Lewis et al. 2002; Gomez et al. 2003; Kauffmann 
et al. 2004), morphology (e.g., Dressier 1980; Goto et al. 
2003; Park et al. 2007), stellar mass (e.g., Kauffmann 
et al. 2004), color (e.g., Hogg et al. 2003; Balogh et al. 
2004; Blanton et al. 2005), and luminosity (e.g., Croton 
et al. 2005; Park et al. 2007) are all strongly correlated 
with local galaxy density. Although there is still debate 
about which, if any, of these relations are “fundamental” 
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(e.g., Hogg et al. 2004; Park et al. 2007), it is clear 
that the mean properties of galaxies we measure depend 
strongly on the type of environment they occupy. 

An obvious first step toward a better understanding of 
environmental processes is to study how they evolve with 
redshift. At higher redshift, the overall population of 
galaxies is younger and have been living within their local 
environment for less time. The environmental processes 
that are most effective and have the shortest timescales 
should be most apparent when comparing galaxies at dif- 
ferent densities in the high redshift universe. The data at 
higher redshift are still somewhat sparse compared to the 
nearby universe but it is beginning to emerge that prop- 
erties such as the SFR (e.g., Elbaz et al. 2007; Cooper 
et al. 2008; Poggianti et al. 2008), color (e.g., Cooper 
et al. 2007) and morphology (e.g., Dressier et al. 1997; 
Postman et al. 2005; Smith et al. 2005, Capak et al. 
2007) are still correlated with local density, albeit differ- 
ently from the nearby universe. 

Of particular interest for understanding environmen- 
tal processes are the cores of rich galaxy clusters. These 
are the most extreme density environments at all red- 
shifts, and if environment is truly an important force in 
galaxy evolution a comparison of the properties of galax- 
ies that live in this environment to those that live in the 
field should provide the largest contrasts. Despite their 
potential value for such studies, and the abundance of 
resources directed at finding distant clusters, there are 
still relatively few confirmed rich clusters at z > 1. 

The major challenge for cluster surveys targeting the 
z > 1 range is the need to be simultaneously deep enough 
to detect either the galaxies or hot X-ray gas in clusters 
and yet wide enough to be able to cover a large area be- 
cause of the rarity of rich clusters at z > 1. The require- 
ment of both depth and area has pushed X-ray detection 
of clusters with current telescopes to the limit. 
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The largest area targeted X-ray cluster surveys are the 
XMM-LSS (Valtchanov et al. 2004; Andreon et al. 2005; 
Pierre et al. 2006) and the XMM-COSMOS (Finoguenov 
et al. 2007), and while these have been successful at dis- 
covering z > 1 clusters (e.g., Bremer et al. 2006), they 
cover areas of only 9 and 2 deg 2 , respectively and are 
therefore limited to fairly low mass systems on average. 
Indeed, X-ray detection of clusters at 2 > 1 is so chal- 
lenging that currently the most promising surveys are 
those searching for clusters serendipitously in the entire 
XMM-Newton archive (e.g., Romer et al. 2001; Mullis 
et al. 2005; Stanford et al. 2006; Lamer et al. 2008). 

Complementary to X-ray detection is optical detec- 
tion of clusters using overdensities of galaxies selected 
using the red-sequence (e.g., Gladders & Yee 2000, 2005; 
Gilbank et al. 2004; Muzzin et al. 2008) or photomet- 
ric redshifts (e.g., Stanford et al. 2005; van Breuklen et 
al. 2007; Eisenhardt et al. 2008). Recently, it has be- 
come clear that the key to discovering clusters above 2 > 
1 with these techniques is the incorporation of infrared 
(IR) data which probes the peak of the stellar emission 
for galaxies at z > 1. 

Although IR surveys have thus far been confined to 
modest areas (ranging from 0.5 - 8.5 deg 2 ) they have 
been extremely successful at detecting z > 1 clusters 
(e.g., Stanford et al. 2005; Brodwin et al. 2006; van 
Breuklen et al. 2007; Zatloukal et al. 2007; Krick et 
al. 2008; Eisenhardt et al. 2008; Muzzin et al. 2008). 
The IR cluster community is now regularly discovering 
clusters at z > 1 and shortly IR-cletected clusters should 
outnumber their X-ray counterparts. 

Currently the largest area IR survey still deep enough 
to detect clusters at z > 1 is the Spitzer Wide-Area 
Infrared Extragalactic Survey (SWIRE, Lonsdale et al. 
2003; Surace et al. 2005). SWIRE covers ~ 50 deg 2 in 
the Spitzer bandpasses and is slightly deeper, and nearly 
a factor of six larger than the next largest IR cluster 
survey, the IRAC Shallow Survey Cluster Search (ISCS, 
Eisenhardt et al. 2008). 

In Muzzin et al. (2008) we demonstrated the poten- 
tial of using the red-sequence method with Spitzer data 
to detect distant clusters using data from the 3.8 deg 2 
Spitzer First Look Survey (FLS, Lacy et al. 2004). In 
2006 we began observations for the Spitzer Adaptation 
of the Red-sequence Cluster Survey (SpARCS), a deep z'- 
band imaging survey of the SWIRE fields. SpARCS aims 
to discover the first large, yet homogeneously-selected 
sample of rich clusters at z >1 using the red-sequence 
method. SpARCS is similar to the RCS surveys (Glad- 
ders & Yee 2005; Yee et al. 2007) which target clusters 
to 2 ~ 1 using an R - z' color except that we use a z' 
- 3.6/im color, which spans the 4000A break at z > 1. 
With a total area effective area of 41.9 deg 2 SpARCS is 
currently the only z > 1 cluster survey that can discover 
a significant number of rare rich clusters. These clusters 
will be extremely valuable for quantifying the evolution 
of galaxy properties in the densest environments at high 
redshift. 

This paper is organized as follows. In § 2 we provide a 
brief overview of the northern component of the SpARCS 
survey (the southern component is summarized in Wil- 
son et al. 2008). In § 3 we discuss the selection of cluster 
candidates that were chosen for followup spectroscopy, 
and in § 4 we present spectroscopic confirmation of two 


z > 1 clusters from early SpARCS data. In § 5 we present 
the dynamical analysis of the clusters followed by a dis- 
cussion of the cluster properties in § 6. We conclude with 
a summary in § 7. 

Throughout this paper we assume an ft m = 0.3, Qa = 
0.7, Ho = 70 km s -1 Mpc -1 cosmology. All magnitudes 
are on the Vega system unless indicated otherwise. 

2. THE SPARCS-NORTH SURVEY 

The SWIRE survey is located in six fields and con- 
tains ~ 50 cleg 2 of imaging in the four IR AC bandpasses 
(3.6pm, 4.5pm, 5. 8pm, and 8. 0pm) and the three MIPS 
bandpasses (24pm, 70pnr, and 160pm). Three of the 
fields are located in the northern hemisphere (ELAIS- 
Nl, ELAIS-N2, and the Lockman Hole), two of the 
fields are located in the southern hemisphere (ELAIS- 
Sl, and the Chandra-S), and one of the fields is equa- 
torial, the XMM-LSS field. A thorough discussion of 
the data reduction, photometry, cluster finding, and the 
SpARCS catalogue for all fields will be presented in a 
future paper by A. Muzzin et al. (2009, in prepara- 
tion) . Here we present a brief summary of the z'-band ob- 
servations of the ELAIS-N1, ELAIS-N2, Lockman Hole, 
and XMM-LSS 14 fields obtained with CFHT /MegaCam; 
hereafter the SpARCS-North Survey. Observations of 
the ELAIS-S1 and Chandra-S fields were obtained with 
the CTIO/MOSAICII and are outlined in the companion 
paper by Wilson et al. (2008). 

The IRAC imaging of the ELAIS-N1, ELAIS-N2, Lock- 
man, and XMM-LSS fields covers areas of 9.8, 4.5, 11.6, 
and 9.4 deg 2 , respectively. In Figure 1 we plot the IRAC 
3.6pm mosaics for these fields. The superposed white 
squares represent the locations of the CFHT /MegaCam 
pointings. The pointings were designed to maximize the 
overlap with the IRAC data, but to minimize the over- 
all number of pointings by omitting regions that have 
little overlap with the IRAC data. There are a total of 
12, 5, 15, and 13 MegaCam pointings in the ELAIS-N1, 
ELAIS-N2, Lockman Hole and XMM-LSS fields, respec- 
tively. 

We obtained observations in the z'-band with 6000s 
of integration time for each pointing in the ELAIS-N1, 
ELAIS-N2 and Lockman Hole fields in queue mode using 
CFHT /MegaCam which is composed of 36 4096 x 2048 
pixel CCDs, and has a field of view (FOV) of ~ 1 deg 2 . 
Omitting the large chip gap areas and regions contami- 
nated by bright stars, the total overlap region with both 
z' and IRAC data is 28.3 deg 2 . 

Photometry was performed on both the z' and 
IRAC mosaics using the SExtractor photometry package 
(Bertin & Arnouts 1996). Colors were determined us- 
ing 3 IRAC pixel (3.66") diameter apertures. The IRAC 
data was corrected for flux lost outside this aperture due 
to the wings of the PSF using aperture corrections mea- 
sured by Lacy et al. (2005). Total magnitudes for the 
IRAC photometry were computed using the method out- 
lined in Lacy et al. (2005) and Muzzin et al. (2008). The 
5(7 depth of the z' data varies depending on the seeing 
and the sky background; however, the mean 5cr depth for 
extended sources is z' ~ 23.7 Vega (24.2 AB). 

3. CLUSTER SELECTION 

14 The XMM-LSS data was obtained as part of the CFHT 
Legacy Survey 
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Fig. 1.- The 3.6/rm mosaics for the four SWIRE fields observable from the northern hemisphere. The location of the SpARCS z'-band 
CFHT/MegaCam pointings are overplotted as white boxes. Each MegaCam pointing covers ~ 1 deg 2 . The observations of the XMM-LSS 
field were obtained as part of the CFHTLS-wide survey. Excluding areas masked by bright stars and missed by MegaCam chip gaps there 
are 28.3 deg 2 with both z' and 3.6/rm observations in the northern fields that can be used for cluster finding. 


Clusters are found in the data using the cluster red- 
sequence algorithm developed by Gladders & Yee (2000; 
2005). Muzzin et al. (2008) used a slightly modified ver- 
sion of the algorithm to detect clusters at 0 < 2 < 1.3 in 
the FLS using an R - 3.6/un color. We use the Muzzin et 
al. (2008) code for the SpARCS data. The change from a 
R - 3.6/irn color to a z' - 3.6/tm is optimum for targeting 
clusters at z > 1, where the z' - 3.6/nn color spans the 
4000A break. Other than the change of using a different 
optical band, the SpARCS algorithm is identical to that 
presented in Muzzin et al. (2008) and we refer to that 
paper for further details of the cluster finding technique. 

After the first semester of z' observations were 
complete there were ~ 14 deg 2 of data with both 
z' and 3.6/im data. From this area we selected 


two rich cluster candidates, both from the ELAIS- 
N2 field with red-sequence photometric redshifts 15 
of z >1.2 for spectroscopic followup. These 
two cluster candidates, SpARCS J163435+402151 
(R.A.: 16:34:35.0, Decl:+40:21:51.0), and SpARCS 

J163852+403843 (R.A.:16:38:52.0, Decl:+40:38:43.0), 
have richnesses, parameterized by Bg Cj #, of 1053 ± 278 
Mpc 18 and 988 ± 270 Mpc 18 , respectively. For a dis- 
cussion of B gc and B gc ^ as cluster richness estimates see 
Yee & Lopez-Cruz (1999), and Gladders & Yee (2005). 
Based on the empirical calibration of B gc vs. M 200 de- 
termined by Muzzin et al. (2007) in the K-band for the 
CNOC1 clusters at z ~ 0.3, this implies M 2 oo = 5.7 x 10 14 

15 Based on a z/ = 2.8 single-burst Bruzual & Chariot (2003) 
model 
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Fig. 2. — Spectra for a subsample of seven galaxies in the clus- 
ter SpARCS J163435+402151. The spectra have been smoothed 
with a 7-pixel boxcar so the sampling matches the instrumental 
resolution. 

Mq and 5.1 x 10 14 M 0 for SpARCS J163435+402151 
and SpARCS J163852+403843, respectively. Although 
Muzzin et al. (2007) found there is a fairly large scatter 
in the B flc vs. M 200 relation, the high richnesses im- 
ply that these candidates are likely to be massive, high- 
redslrift systems. 

4. SPECTROSCOPIC DATA 

Multislit nod-and-shuffle (N&S) spectroscopy of 
galaxies in SpARCS J163435+402151, and SpARCS 
J163852+403843 were obtained using GMOS-N on Gem- 
ini as part of the program GN-2007A-Q-17. We used 
the R150 grating blazed at 7170A with 1" width slits. 
This provided a resolving power of R = 631 which corre- 
sponds to a resolution of ~ llA, or ~ 250 km s _1 at the 
estimated redshift of the clusters. For all observations 
we used 3" long microslits, corresponding to roughly 
four seeing-disks, allowing a two seeing-disk spacing be- 
tween the nod positions. We observed three masks for 
SpARCS J163435+402151 and four masks for SpARCS 
J163852+403843. One mask for each of the clusters was 
observed in “micro-shuffle” mode, but the majority were 
observed in “band-shuffle” mode. All masks were ob- 
served using the RG615 filter which blocks light blueward 
of 6150A so that multiple tiers of slits could be used. 

Unlike micro-shuffle where the shuffled charge is stored 
directly adjacent to the slit location, band-shuffle shuffles 
the charge to the top and bottom third of the chips for 
storage. While technically it is the least-efficient N&S 
mode in terms of usable area for observations (only the 
central 1.7' of the total 5' FOV can be used) it is ex- 
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Fig. 3. — Same as Figure 2 but for galaxies in the cluster SpARCS 
J 163852+403843. 

tremely efficient for observations of high-redshift clusters 
because it allows the microslits to be packed directly be- 
side each other in the cluster core with no requirement for 
additional space for storing the shuffled charge. In band- 
shuffle mode we were typically able to locate between 
20-26 slits, including three alignment stars, per mask in 
the central 1.7' around the cluster. At z ~ 1.2 the 1.7' 
FOV corresponds to a diameter of 850 kpc, roughly the 
projected size of a massive cluster. 

Slits were placed on galaxies with priorities in the fol- 
lowing order: Priority 1, galaxies with colors ± 0.6 mag 
from the red-sequence and 3.6/rm < 16.9. Priority 2, 
galaxies with colors ± 0.6 mag from the red-sequence, 
3.6/xm > 16.9 and z' < 23.5. Priority 3, galaxies with 
colors > 0.6 bluer than the red-sequence, but <1.0 mag 
bluer and 3.6/im > 16.9 and z' < 23.5. Priority 4, same 
as priority 3 but for galaxies with colors bluer than the 
red-sequence by 1.0-1. 4 mag. Priority 5: all galaxies 
with 23.5 < z' < 24.5. Roughly speaking, Priorities 1 
through 4 can be described as bright red-sequence, faint 
red-sequence, blue cloud, and extreme blue cloud galax- 
ies, respectively. 

For each mask we obtained a total of 3 hrs of inte- 
gration time by combining six exposures with 30 mins 
of integration time. The six frames were obtained using 
15 nod cycles of 60s integration time per cycle. Each of 
the 6 exposures was offset by a few arcseconds using the 
on-chip dithering option. 

4.1. Data Reduction 

Data were reduced using the GMOS IRAF package. 
We subtracted a bias and N&S dark from each frame. 
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The N&S darks are taken using the same exposure times 
and using the same charge shuffling routine as the sci- 
ence observations, but with the shutter closed. Regions 
with poor charge transfer efficiency cause electrons to be- 
come trapped during the repeated charge shuffling used 
in the observations. Such charge traps can be identified 
and corrected using dark frames taken with the same 
N&S settings. Images were registered using bright sky 
lines and sky subtracted using the complementary stor- 
age area using the “gnsskysub” task. Final mosaics are 
made by coadding the sky subtracted images. 

One dimensional spectra were extracted using the 
iGDDS software (Abraham et al. 2004). Wavelength 
calibration for each extracted spectrum was performed 
using bright sky lines from the unsubtracted image, also 
with the iGDDS software. Wavelength solutions typi- 
cally have an rms < 0.5 A. We determined a relative flux 
calibration curve using a long slit observation of the stan- 
dard star EG131. 

Redshifts were determined interactively for each spec- 
trum by comparing with the templates available in 
iGDDS. Most of the redshifts were identified using the 
early-, intermediate-, and late-type composite spectra 
from the Gemini Deep Deep Survey (Abraham et al. 
2004). The final redshifts was determined using the av- 
erage redshift from all absorption and emission lines that 
were detected. The vast majority of redshifts were de- 
termined by identifying the [Oil] 3727A-doublet emission 
line (which is not resolved at our resolution), or the Cal- 
cium II H+K absorption lines. Many of the spectra also 
show the Balmer series lines. We list the spectroscopic 
members of SpARCS J163435+402151 and SpARCS 
J163852+403843 in Tables 1 and 2, and the spectroscop- 
ically confirmed foreground/background galaxies in Ta- 
bles 3 and 4. We also plot examples of some cluster 
galaxy spectra in Figures 2 and 3. R, z', and 3.6 /cm 
color composites of the two clusters are shown in Figures 
4 and 5. The white squares denote the spectroscopically 
confirmed cluster members and the green squares denote 
the spectroscopically confirmed foreground/background 
galaxies. 


5. CLUSTER VELOCITY DISPERSIONS 

For both clusters there are a sufficient number of red- 
shifts to determine a velocity dispersion ( a v ) and there- 
fore a dynamical mass. Our o+’s are determined us- 
ing the method detailed in Blindert (2006). Briefly, we 
make a rejection of near-field non-cluster members using 
a modified version of the Fadda et al. (1996) shifting- 
gap procedure. This method uses both the position and 
velocity of galaxies to reject interlopers. In Figure 6 we 
plot the relative velocities of the cluster galaxies as a 
function of projected radius. Two galaxies in SpARCS 
J163435+402151 are rejected as near-field interlopers 
and three are rejected in SpARCS J163852+403843. The 
rejected galaxies are plotted as crosses in Figure 6 and 
are not used in computing the mean redshift of the clus- 
ter or a v . Once outliers are rejected the redshift of 
the clusters is determined using the remaining galax- 
ies. The spectroscopic redshift of the clusters is 1.1798 
and 1.1963 for SpARCS J163435+40215 and SpARCS 
J163852+403843, respectively. 


TABLE 1 

Spectroscopic Cluster. Members in SpARCS 
J163435+402151 


ID 

R.A. 

J2000 (Deg.) 

Decl. 

J2000 (Deg.) 

z' 

Mag Vega 

Zspec 

Mask 1 

3 

248.6589 

40.35303 

21.54 

1.179 

5 

248.6497 

40.36148 

21.58 

1.181 

6 

248.6467 

40.36418 

20.95 

1.166 

21 

248.6513 

40.35238 

22.14 

1.174 

22 

248.6605 

40.35474 

21.91 

1.178 

25 

248.6510 

40.35855 

21.44 

1.180 

28 

248.6605 

40.36132 

22.19 

1.182 

34 

248.6527 

40.36477 

22.78 

1.182 

37 

248.6496 

40.36696 

21.82 

1.180 

41 

248.6459 

40.36841 

21.97 

1.185 

45 

248.6440 

40.37413 

22.14 

1.187 

3027 

248.6546 

40.35461 

21.59 

1.181 

Mask 2 

20 

248.5991 

40.34870 

23.10 

1.170 

27 

248.6548 

40.35104 

22.23 

1.164 

36 

248.6171 

40.36649 

22.36 

1.184 

57 

248.6700 

40.39159 

22.80 

1.176 

Mask 3 

2 

248.6584 

40.34934 

20.50 

1.175 

38 

248.6474 

40.36152 

21.67 

1.178 

69 

248.6669 

40.39132 

21.89 

1.178 


TABLE 2 

Spectroscopic Cluster Members in SpARCS 
J 163852+403843 


ID 

R.A. 

Decl. 

z' 

Zspec 


.12000 (Deg.) 

.12000 (Deg.) 

Mag Vega 

Mask 1 

5 

249.7132 

40.63952 

21.51 

1.194 

7 

249.7144 

40.64159 

21.68 

1.194 

8 

249.7152 

40.64527 

21.41 

1.195 

1007 

249.6974 

40.63383 

22.50 

1.199 

1010 

249.7151 

40.63804 

22.48 

1.202 

1016 

249.7006 

40.64276 

22.20 

1.200 

1025 

249.7137 

40.64952 

22.46 

1.200 

1028 

249.7104 

40.65532 

21.96 

1.190 

1031 

249.7029 

40.65968 

22.85 

1.176 

3026 

249.6668 

40.64420 

22.73 

1.195 

Mask 2 

1018 

249.7032 

40.64533 

22.67 

1.186 

1020 

249.7099 

40.64631 

22.39 

1.202 

1024 

249.7135 

40.64828 

22.24 

1.198 

1026 

249.7126 

40.65311 

22.26 

1.200 

1030 

249.7581 

40.63628 

21.95 

1.195 

2016 

249.6610 

40.65976 

22.22 

1.197 

Mask 3 

1017 

249.7136 

40.64511 

22.38 

1.198 

1019 

249.7163 

40.64602 

22.29 

1.198 

2009 

249.7236 

40.63171 

22.10 

1.188 

4029 

249.7619 

40.65497 

23.67 

1.195 

Mask 4 

2 

249.6947 

40.61541 

21.24 

1.196 

3 

249.6986 

40.62432 

21.73 

1.200 

10 

249.7314 

40.66114 

21.25 

1.192 

11 

249.7378 

40.66462 

21.72 

1.194 

13 

249.6680 

40.67089 

21.17 

1.194 

2014 

249.7532 

40.64920 

22.58 

1.196 

2030 

249.7049 

40.68847 

21.87 

1.175 

1012 

249.6992 

40.63979 

23.06 

1.196 

1021 

249.7217 

40.64632 

21.90 

1.195 

3019 

249.7189 

40.63320 

23.17 

1.219 

3046 

249.7509 

40.68401 

22.05 

1.172 
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Fig. 4. — Left: Rz , 3.6//rn color composite of the cluster SpARCS J163435+402151 at z = 1.1798. The R and z' images have been 
convolved to match the 3.6pm PSF. The FOV of the image is ~ 3.5' across. Right: Same as left panel but with spectroscopically confirmed 
cluster members marked as white squares and spectroscopically confirmed foreground/background galaxies marked as green squares. 




Fig. 5. — As Figure 4, but for the cluster SpARCS J163852+403843 at z = 1.1963. The FOV of the images is ~ 4.5 / across. 


TABLE 3 

Spectroscopic Foreground/Background Galaxies 
in Field of SpARCS J163435+402151 


TABLE 4 

Spectroscopic Foreground/Background Galaxies 
in Field of SpARCS J163852+403843 


ID 

R.A. 

Deck 

t! 

Zspec 


J2000 (Deg.) 

J2000 (Deg.) 

Mag Vega 

Mask 1 

4 

248.6129 

40.35610 

20.89 

1.108 

26 

248.6600 

40.35898 

21.38 

1.004 

31 

248.6042 

40.36213 

21.66 

1.105 

45 

248.6823 

40.38687 

23.64 

0.925 

Mask 2 

23 

248.60280 

40.33815 

23.79 

1.337 

24 

248.60710 

40.33762 

24.05 

1.255 

Mask 3 

i 

248.67280 

40.33252 

21.16 

1.348 

15 

248.65750 

40.33183 

21.32 

0.780 

65 

248.67920 

40.38581 

23.45 

1.108 

67 

248.63340 

40.38773 

23.05 

0.811 


The cr v ’ s are determined using the “robust” estima- 
tor suggested by Beers et al. (1990) and Girardi et al. 
(1993). The robust estimator is simply the biweight es- 


ID 

R.A. 

Decl. 

z ' 

Zspec 


J2000 (Deg.) 

J2000 (Deg.) 

Mag Vega 

Mask 1 

3021 

249.7659 

40.63502 

22.14 

0.670 

Mask 2 

3024 

249.7503 

40.64119 

23.61 

0.875 

4020 

249.7305 

40.63887 

23.98 

1.391 

4030 

249.7666 

40.65555 

23.56 

0.776 

Mask 3 

3034 

249.6953 

40.65273 

21.96 

1.393 

3036 

249.6603 

40.65813 

23.47 

1.386 

Mask 4 

1 

249.7320 

40.60904 

22.10 

0.963 

12 

249.6992 

40.66802 

20.96 

0.784 

2003 

249.6765 

40.60412 

22.44 

1.017 

2024 

249.7496 

40.67251 

23.50 

0.848 

2026 

249.6744 

40.67759 

21.46 

0.771 

3004 

249.7519 

40.60677 

21.44 

0.768 

3041 

249.7611 

40.67033 

21.17 

0.784 
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timator for systems with > 15 members, and the gap- 
per estimator for systems with < 15 members. As dis- 
cussed in those papers and Blindert (2006) these estima- 
tors are more robust than standard deviations as they 
are less sensitive to outliers, which may still persist even 
after the initial shifting-gap rejection. Using the “ro- 
bust” estimator, SpARCS J163435+402151 and SpARCS 
J163852+403843 have u v = 490 ± 140 km s" 1 and 650 
± 160 km s -1 , respectively, where the errors have been 
determined using Jackknife resampling of the data. 

We estimate the dynamical mass using M 2 oo, the mass 
contained within r 2 oo, the radius at which the mean in- 
terior density is 200 times the critical density (p c ). We 
use the equation, 

A/200 = ' 200p c , (1) 

with the dynamical estimate of r 2 oo from Carlberg et al. 
(1997), 

_ >/3 a 

r20 ° - 1 may (2) 

where H(z) is the Hubble constant at the redshift of 
the cluster. From these relation we derive r 2 oo = 
0.62 ± 0.18 Mpc, and 0.82 ± 0.20 Mpc for SpARCS 
J163435+402151 and SpARCS J163852+403843, respec- 
tively. From Equation 2, these imply M 2 oo = (1-0 ± 
0.9) x 10 14 M 0 and (2.4 ± 1.8) x 10 14 M 0 for SpARCS 
J163435+402151 and SpARCS J163852+403843, respec- 
tively. 

6. DISCUSSION 

6.1. Red-sequence Photometric Redshifts 

In Figure 7 we plot the z! - 3.6/rm vs. 3.6/xm color mag- 
nitude relation for galaxies at projected radii (R) < 550 
kpc in the fields of both clusters. Spectroscopically con- 
firmed members and confirmed foreground/background 
galaxies are plotted as red and blue diamonds, respec- 
tively. The dotted line in both panels is the best fit 
line with the slope fixed at zero, to the spectroscop- 
ically confirmed members. These lines indicate that 
the red-sequence galaxies have z' - 3.6/zm colors of 4.77 
and 4.82 for SpARCS J163435+402151 and SpARCS 
J163852+403843, respectively. Using a solar metallic- 
ity, Bruzual & Chariot (2003) simple stellar population 
(SSP) with a Zf = 4.0 these colors imply photometric 
redshifts of 1.19 and 1.20, in excellent agreement with 
the spectroscopic redshifts. At z ~ 1.2, the red-sequence 
photometric redshifts do not depend strongly on the cho- 
sen Zf. If we instead use a z/ = 2.8 SSP, the red-sequence 
color would predict redshifts of 1.21 and 1.24, and for a 
Zf = 10.0 SSP it would predict redshifts of 1.13 and 
1.15. However, the color differences between all these 
models at fixed redshift are small (< 0.1 mag), and so it 
is not possible to distinguish between different formation 
epochs without more data. Still, the close agreement 
between the red-sequence photometric redshift derived 
using a reasonable Zf and the spectroscopic redshift is 
encouraging for the use of red-sequence photometric red- 
shifts for clusters at z > 1. 

6.2. Mass vs. Richness 

Both clusters have lower masses than predicted by 
their richness by factors of ~ 6 and 2 for SpARCS 


J163435+402151 and SpARCS J163852+403843, respec- 
tively, although due to the large error bars the differ- 
ences are only significant at ~ 1 and 2a, respectively. 
Whether this represents a redshift evolution in the B sc - 
M 20 o scaling relation, or is simply a richness-selected Ed- 
dington bias 16 is impossible to determine using only two 
clusters. Both Gilbank et al. (2007) and Andreon et 
al. (2007) found that for a small sample of clusters at 
z ~ 1 the cluster richnesses were still consistent with 
their velocity dispersions based on relations calibrated 
at lower redshift, although both parameters have large 
uncertainties in their measurements. More clusters with 
well-determined a v and P> gc will be needed to test if the 
cluster scaling relations at z > 1 are similar to those at 
lower redshift. 

7. SUMMARY 

We have presented a brief summary of observations for 
the northern component of the SpARCS survey. Using 
Gemini N&S spectroscopy we confirmed two rich cluster 
candidates at z ~ 1.2 selected from early survey data. 
We find that the photometric redshifts from the color 
of the cluster red-sequence agree extremely well with 
the spectroscopic redshifts. Both clusters have a smaller 
M 20 o than would be expected from their richness if we 
use the B ffc - M 20 o scaling relation calibrated at z ~ 0.3. 
Whether this represents a true evolution in the cluster 
scaling relations at z > 1.2 or is simply a selection bias 
will require well-determined M 2 oo for a larger sample of 
clusters. 

Overall, the confirmation of both SpARCS 
J163435+402151 and SpARCS J163852+403843 as 
bona fide massive clusters at z > 1 provide strong 
evidence that the red-sequence technique is an effective 
and efficient method for detecting clusters at z > 1 
(see also Wilson et al. 2008 who present a confirmed 
z = 1.34 cluster from the southern component of the 
SpARCS survey). The complete SpARCS catalogue 
contains hundreds of cluster candidates at z > 1 and 
promises to be one of the premier data sets for the study 
of cluster galaxy evolution at z > 1. 

Based on observations obtained with 
MegaPrime/MegaCam, a joint project of CFHT 
and CEA/DAPNIA, at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National 
Research Council (NRC) of Canada, the Institut Na- 
tional des Sciences de l’Univers of the Centre National 
de la Recherche Scientifique (CNRS) of France, and the 
University of Hawaii. This work is based in part on data 
products produced at TERAPIX and the Canadian 
Astronomy Data Centre as part of the Canada-France- 
Hawaii Telescope Legacy Survey, a collaborative project 
of NRC and CNRS. 

Based on observations obtained at the Gemini Obser- 
vatory, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under a cooper- 
ative agreement with the NSF on behalf of the Gemini 

16 We followed up two of the richest clusters in our early dataset. 
The cluster mass function is steep at high redshift and low mass 
systems greatly outnumber high mass systems. Due to scatter 
in the mass-richness relation lower mass systems with abnormally 
high richnesses may be more common than truly massive systems. 
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Fig. 6. — Left Panel: Galaxy velocities relative to the cluster mean velocity as a function of radius for SpARCS J163435+402151. Right 
Panel: Same as left panel but for SpARCS J163852+403843. Galaxies marked with an “x” are more likely to be near-field objects than 
members of the cluster and are not used in the computation of the velocity dispersion. 




Fig. 7. — Left Panel: z' - 3.6pm vs. 3.6pm color magnitude diagram for galaxies at R < 550 kpc in the field of the cluster SpARCS 
J163435+402151. Spectroscopically confirmed cluster members and foreground/background galaxies are plotted as red and blue diamonds, 
respectively. The dotted line is the best fit line to the confirmed cluster members with the slope fixed at zero. Right Panel: Same as left 
panel but for SpARCS J 163852+403843. 


partnership: the National Science Foundation (United 
States), the Science and Technology Facilities Coun- 
cil (United Kingdom), the National Research Council 
(Canada), CONICYT (Chile), the Australian Research 
Council (Australia), Ministerio da Ciencia e Tecnologia 


(Brazil) and SECYT (Argentina) 

This work is based in part on observations made with 
the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology under a contract with NASA. 
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